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Abstract. The rotational spectra of DCO''' and its two isotopomers D^'^CO^and 
DC^**0 , produced in a negative glow discharge cell, have been recorded in the 
137-792 GHz region, which includes lines from J = 2 ^ 1 up to J = 11 ^ 10. 
The determined rotational and centrifugal distortion constants allow to predict 
the DGO"'' rotational spectrum up to 1000 GHz with an accuracy of 1 part in 10^ 
or better. This is important for kinematic studies of dense molecular cloud cores 
and for future far-infrared observations. We also report on the first detection of 
the hyperfine structure of the DCO^(l-O) line made at the IRAM 30m antenna, 
toward the quiescent starless cloud core L1512, in the Taurus Molecular Cloud. 
We point out that this is the first observation of the hyperfine splitting due to the 
deuteron. This allowed us to quantify the effects of the hyperfine splitting on the 
line width determination; if the hyperfine structure is not taken into account in 
the line fit, the DCO^(l-O) line width is overestimated by a significant factor 
2). 
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molecules - radio lines: ISM 

1. Introduction 



+ 



In the past few years it has been realized how important is to know with high preci- 
sion the frequencies of the molecular transitions used to investigate the physics and the 
chemistry of interstellar clouds and star forming regions (e.g. Mardones et al. 1997; Lee, 
Myers & Tafalla 1999, 2001). This is because spectral line observations are unique tools 
to study internal motions of dense molecular cloud material, within which stars will be or 
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are forming. In this sense, molecular line observations are crucial to test current models 
of cloud core and star formation and give insights on molecular cloud evolution. 

In particular, before low-mass stars form, the progenitor (or starless) core is typically 
a cold (~ 10 K) and quiescent region, where the (optically thin) line widths are slightly 
(approximately a factor of 2) larger than the thermal value Ai^t (e.g. Caselli et al. 2002). 
This means that Auobs ~ 2 x Avt ^ 0.3 (r/lOK)!/^ (v4x/30 amu)"i/2 knis-\ where 
T is the gas temperature and Ax is the molecular weight. This already shows that to 
determine the LSR velocity of a cold (T = 10 K) cloud with an accuracy comparable to the 



standard deviation of the corresponding gaussian line profile (ctvlsr ~ Avo^b/ \/Sln{2) 5s 
0.1 kms^^), for a molecular species with Ax ~ 30 amu (such as CO, HCO^, and DCO^), 
the observed line frequency should have an uncertainty of Aj/(kHz) ^ 0.1 x 3.3 x i'o(GHz) 
(where z/q is the rest frequency), or about 30 kHz at millimeter wavelengths. 

If observed lines are optically thick and if they are tracing inward motions, they 
appear broader and asymmetric, typically double peaked, with the blue peak stronger 
than the red peak (e.g. Leung & Brown 1977; Zhou et al. 1990). However, to quantify 
the infall velocity one needs to observe an optically thin tracer and compare the thin and 
thick profiles (e.g. Myers et al. 1996). For this purpose, one typically uses the normalized 
velocity difference between the thin and thick lines {6V, as defined by Mardones et al. 
1997). In this case, the condition SV/agv > 3 is only obtained if the error on V is cry 
SsO.Ol kms~^ (assuming average values of dV ^ 0.2, aAv ^ 0.03 kms~^, and AV ~ 
0.3 kms~^; see Mardones et al. 1997). This implies a frequency precision of (Ty(kHz) ^ 
0.03xt'o(GHz), or about 3 kHz at ^^100 GHz. All this demonstrates the need of laboratory 



work to accurately measure line frequencies. Gottlieb et al. (2003) have also pointed out 
that it is a challenge to the spectroscopist the determination of line frequencies with a 
precision large enough to identify motions in molecular cloud cores from tiny velocity 
shifts, which typically require a precision of a few parts in 10^. In the case of DCO^ this 
is equivalent to 1.7% of the Doppler line width at 77 K. 

DCO^ is a particularly suitable species to study the chemistry, in particular the 
deuterium fractionation (e.g. Saito et al. 2002) and the electron fraction (e.g. Guelin et 
al. 1977, Caselli et al. 1998, Williams et al. 1998) and the kinematics of the dense cores. 
The lower rotational transitions are easily detected, as it is proved by the large amount 
of literature on this molecular ion, from interstellar clouds (e.g. Wootten et al. 1982) to 
dense cloud cores (e.g. Butner et al. 1995), to circumstellar disks (van Dishoeck et al. 
2003). DCO^ has been recently surveyed in a large sample of starless cores by Lee et al. 
(2004), with the aim of detecting infall motions. However, the previous measurements of 
DCO^millimeter-wave frequencies, carried out in two different laboratories (Bogey et al. 
119811 ISastry et al. 1981| ), gave results differing up to 45 part in 10® for the J = 3 ^ 2 
line (see Table therefore there is a need for more accurate frequency measurements of 
DCO^ lines in order to use them in the study of star-forming molecular cloud cores. In 
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this paper, we repeat with far better precision the old measurements of milHmeter-wave 
transitions of DCO+, D^^CO^, and DC^^^O^ (Bogey at al. lTMUl . and we add new mea- 
surements of submiUimeter-wave transitions: the spectroscopic constants derived from 
the analysis of the frequency data enable to predict spectra up to 1000 GHz with an 
accuracy of at least 1 part in 10* for DCO^, 5 parts in 10* for D^^CO^, and of 1 part 
in 10^ for DC^*0^; all the isotopomers have been observed in natural abundance. 



Table 1. Observed rotational transition of DCO^(MHz) 



J' 


J 


Boeev at al.Il^STl |Sastry at al. 1981| 


diff.^/ppb 


this work 


radiastronomy*^ 


2 


1 


144 077.342 


144 077.319 


160 


144 077.2890 


144 077.321 


3 


2 


216112.701 


216112.604 


449 


216112.5822 




4 


3 


288 144.018 


288 143.911 


371 


288 143.8583 


288 143.911 


5 


4 


360 169.980 


360 169.881 


275 


360 169.7783 


360 169.830 


6 


5 




432 189.033 




432 189.0052 





" difference between measurements of Bogey at al. 119811 and of |Sastry et al. 1981| 
ILovas 20041 



2. Observations 

2.1. Laboratory measurements 

The laboratory spectrum was observed with a frequency-modulated millimeter-wave 
spectrometer l|Cazzoli fc Dore 1990a|l equipped with a double-pass negative glow dis- 
charge cell IjDore et al. 1999|) made of a Pyrex tube 3.25 m long and 5 cm in diame- 
ter. The radiation source was a frequency multiplier driven by a Gunn diode oscilla- 
tor working in the region 68-76 GHz (Farran Technology Limited) for the lines up to 
J = 4 <~ 3; for the high frequency transitions, a doubler in cascade with a multiplier 
(RPG - Radiometer Physics GmbH) was driven by Gunn oscillators working in the re- 
gion 75-115 GHz (J. E. Carlstrom Co and RPG) to cover frequencies up to 792 GHz. 
Two phase-lock loops allow the stabilization of the Gunn oscillator with respect to a 
frequency synthesizer, which is driven by a 5-MHz rubidium frequency standard. The 
frequency modulation of the radiation is obtained by sine- wave modulating with low dis- 
tortion (total harmonic distortion less than 0.01%) the reference signal of the wide-band 
Gunn-synchronizcr. The signal, detected by a liquid-hclium-cooled InSb hot electron 
bolometer (QMC Instr. Ltd. type QFI/2), is demodulated at 2-/ by a lock-in amplifier. 

DCO"'' was produced by flowing a 1:1 mixture of CO2 and D2 (1 mTorr) in Ar buffer 
gas with a total pressure of about 5 mTorr and discharging with a current of a few mA; 
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D^'^CO^and DC^^O^were observed in natural abundance. The cell was cooled at 77 K 
by liquid nitrogen circulation, and an axial magnetic field up to about 400 G was applied 
throughout the length of the discharge. With this longitudinal magnetic field applied, 
ions are produced and observed in the negative glow IjDe Lucia et al. 1982|l . which is a 
nearly field free region and where they are expected to show negligible Doppler shift 
due to the drift velocity, occurring, instead, in the positive column ( |Sastry et al. 1981) ) 
where a low axial electric field is present. In addition, the double-pass arrangement would 
compensate for such a shift, if present. 

A typical spectrum is recorded by sweeping the frequency up and down (several 
times if signal averaging is needed) in steps of 5 kHz at a rate of 0.8 MHz s~^, with 
a lock-in amplifier time constant of 3 ms and a frequency modulation depth compara- 
ble to the Doppler width or larger for the rarer isotopomers. Since we have full flex- 
ibility in controlling scanning rate, number of data points and modulation depth, the 
values of these parameters have been adjusted to prevent any bias of the measured tran- 
sition frequency, which is recovered from a line shape analysis of the spectral profile 
HCazzoh fc Dore IQQOairPore 2003|l . 

2.2. DCO^ (1-0) toward the quiescent cloud core L1512 



The DCO^(l-O) spectrum toward the quiescent Taurus starless core L1512 (Fig. ^ 
has been obtained with the IRAM-30m antenna, located at Pico Veleta (Spain) in August 
2004. The adopted coordinates were R.A.(1950) = 05*^ 00'" 54.4% Dec.(1950) = 32° 39' 
00.0". We used the recent extension of the 3mm tuning range below 80 GHz and carried 
out the observations using the frequency switching technique with a throw of 3.9 MHz. 
The spectral resolution is 3.3 kHz, which corresponds to 14 m/s at 72 GHz. Given that 
below ~ 75 GHz the image gain depends strongly on the tuning parameters and often 
varies across the bandpass of the receiver, we measured the sideband gain ratio and 
corrected the observed spectrum accordingly. The pointing was checked every 2 hours 
and found to be accurate within about 10" because of anomalous refraction problems 
during very good atmospheric conditions. The half power beam width (HPBW) was 33". 
The units are main beam brightness temperature, assuming a source filling factor of 
unity. The beam and forward efficiencies were 0.79 and 0.95, respectively. The rms noise 
of the spectrum is 0.09 K. 

Figure n clearly shows that the DGO^(l-O) line profile is complex and contains struc- 
ture. Indeed, the three detected features are consistent with the hyperfine splitting due to 
the interaction between the molecular electric field gradient and the electric quadrupole 
moments of the deuteron (spin =1), which produces the line splitting into three compo- 
nents. The velocity shifts predicted by theory are larger than the observed separations 
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between the three components by about 0.1 kms^^, suggestmg that the quadrupole con- 
stant needs to be refined. Thus, we estimated the separations performing a gaussian fit to 
the three components of the DCO^(l-O) spectrum. This allowed us to derive the velocity 
separations of -0.198±0.012 kms^^ and 0.256±0.017 kms^^ between the main (F',F = 
2,1) component and the low (F',F = 1,1) and high (F',F = 0,1) velocity components, 
respectively. With this new values of the frequency shifts and the relative intensities of 
the three components (3/9,5/9,1/9 for the (1,1), (2,1) and (0,1) hyperfines, respectively) 
we then performed an hfs fit using the IRAM reduction package CLASS, shown in Figure 
1 (upper panel) by the thick curve. Forcing the resulting LSR velocity (Vlsr) to coincide 
with that measured with the high sensitivity N2H+(l-0) spectrum reported in Caselli et 
al. (1995), which has been corrected to account for the new N2H+(l-0) frequency value 
(93176.2608 MHz ± 6 kHz for the component Fi,F = 0,1^1,2) estimated by Dore et al. 
(2004), we obtain the frequency values reported in Table 2. These new frequencies have 
finally been used to refine the quadrupolar constant (see next section). 

The hfs fit in CLASS gives the following parameters: Vlsr = 7.094±0.004 kms"\ 
intrinsic (i.e. corrected for optical depth effects) line width Av — 0.149±0.008 kms^^, 
total optical depth (i.e. the sum of the optical depths of the three hyperfines) ttot 
= 1.9±0.5, excitation temperature T^x — 4±1 K. We note here that the observed line 
width is only 1.25 times larger than the DCO^ thermal width Avt, assuming a kinetic 
temperature of 10 K (or Av — Avt if T = 15 K). If the hyperfine structure is not taken 
into account, the DCO^(l-O) line width obtained with a simple Gaussian fit is a factor 
of 2.3 times larger, in net contrast with the narrow (~ 0.18 kms^^) N2H+(l-0) line 
width observed in the same position (Caselli et al. 1995). Therefore, we conclude that 
the differences in the widths of DCO^(l-O) and NH3(1,1) lines observed in previous work 
(e.g. compare Butner et al. 1995 with Benson & Myers 1989) are likely due to neglecting 
the hyperfine splitting. The bottom panel of Fig. ^ shows how the line profile changes 
in case of large optical depths. For ttot = 10, the weakest (F', F = 0,1) component 
becomes bright enough to further enlarge the line width by other ^0.25 kms^^ (so that 
the total linewidth will be about 0.6 kms^^, if the splitting is not considered in the fit). 

The effects of the hyperfine structure on the line width and profile become less impor- 
tant for higher J transitions. In Fig. El the simulated spectrum of the DC0^(2-1) line, 
with the same parameters as the DCO'''(1-0) line in Fig. ^ is shown. The six hyperfine 
components are blended together, but there is still some (^^20 %) extra-broadening and 
a small (0.01 km s^^) line center shift resulting from the hfs splitting, which should be 
taken into account when analysing DC0^(2-1) profiles. In particular, this may be im- 
portant to quantify the infall velocity in the central regions of starless cores, as recently 
done by Lee et al. (2004). On the other hand, the components of the J = 3^2 transition 
(as well as those of higher J transitions) are heavily blended together, and the line width 
measurement is not affected at all by their presence. 
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Vlsr (km/s) 

Fig. 1. ( Top panel) Observed spectrum of DCO"*" (1-0) toward the quiescent starless cloud 
core L1512, in the Taurus Molecular Cloud. The units are in main beam brightness tem- 
perature, assuming a unity filling factor. This is the first spectrum showing the hyperfine 
structure of the line due to the deuteron electric quadrupole interaction. The line is sig- 
nificantly broadened by the hyperfine structure, and the width is a factor of 2.3 larger if 
the splitting is not taken into account in the fit. The black curve is the hfs fit performed 
in CLASS. {Bottom panel) Simulated spectrum of the DCO'''(1-0) line assuming the 
same parameters of the observed line (Tex = 4 K and intrinsic At; = 0.15 kms~^) but 
with different values of ttot (from least to most intense curves: 0.1, 1, 2 (as observed), 
5, and 10. Vertical bars indicate the positions of the three hyperfines with their relative 
intensities in the optically thin limit. 
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6.5 



7 

Vlsr (km/s) 



7.5 



Fig. 2. Simulated spectrum of the J — 2— >1 transition of DCO in L1512, including the 
hyperfine splitting due to the deuteron electric quadrupole interaction (black curves). 
The assumed excitation temperature is 4 K and the intrinsic linewidth is 0.15 km s^^, 
as found for the DCO^(l-O) line, whereas the total optical depth is 1 (thick curve) and 
10 (thin curve). Dashed curves are the expected DC0^(2-1) lines in case of no hyperfine 
structure, again assuming ttot ~ 1 (thick dashed curve) and 10 (thin dashed curve). 
Vertical lines mark the velocities of the 6 hyperfine components with relative intensities 
given by their statistical weight. The hyperfine splitting leads to a '^20 and ~30% increase 
of the line width for ttot = 1 and 10, respectively. 

3. Analysis and discussion 

The measured transition frequencies of DCO^, D^'^CO^, and DC^^O^are listed in 
Tables El El and El they are mean values obtained from 5 to 14 measurements. The 
standard errors of the mean result unrealistically small (for DCO^they range from 0.1 
to 1 kHz), therefore the uncertainties reported in the last column of each table were 
estimated from the ratio of linewidth to signal-to-noise ratio, and were roughly scaled 
according to the standard errors. 

The experimental data were fitted, in a weighted least-squares procedure, to the 
standard expression for the frequency of the rotational transition J + 1 ^ J: 



1^0 = 2B{J + 1) - 4Lij(J + 1)=^ -t- Hj{J + if [{J + 1)=* - J\ 



(1) 
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where B is the rotational constant and Dj and Hj arc the quartic and sextic centrifugal 
distortion constants, respectively; the weights were the inverse-square of the uncertainties. 
The standard deviation (jfn of the fit for each isotopomer is just a few kHz with the 
inclusion of Hj as fit parameter. An assessment of the significance of the latter can 
be carried out by comparison with DCN centrifugal distortion constants: the values 
IjBriinken et al. 2004|l reported in Tables [3 and |21 result comparable with the values of 
Dj and Hj derived in this work for DCO^'^and D^^CO^, therefore these constants may 
be assumed correctly determined. 

In the case of DCO^, the three hyperfine frequencies of the J — I ^ transition, 
accurately determined as reported above, could be used to derive the quadrupole coupling 
(eqQ) and spin rotation (C/) constants of the D nucleus (spin quantum number / = 1) 
according to the expression: 

i^hf^'yQ- eqQ[Y{J + 1, /, F') - Y{J, I, F)] + Ci[C[J + 1, J, F') - C{J, J, F)], (2) 

where F is the total angular momentum [F = J + I) quantum number, C{J,I,F) — 
F{F + 1) - /(/ + 1) - J{J + 1) and the Casimir function Y{J, I, F) is given by: 

^ ' ' 2(2/- 1)(2J + 3)7(2/ - 1) ^ ' 

The Pickett's SPFIT fitting program fPickett allowed us to carry out simultane- 

ously the hfs and centrifugal analyses in a global fit of laboratory and radioastronomical 
transition frequencies, whose weights were the inverse-square of the uncertainties. For 
the hyperfines, their assumed uncertainties account for both the standard errors of the 
hfs fit and the 6 kHz uncertainty of the N2H+(l-0) reference line IjDore et al. 2004|l . 

Tables 121 13 and 01 list, for each isotopomer, the experimental frequencies with their 
uncertainties and the residuals of the least-squares fit carried out to determine the spec- 
troscopic constants reported at the bottom of the table. These accurate constants (see 
next paragraphs for a discussion of accuracy) allow us to recommend the rest frequen- 
cies reported in Table the hyperfine components of the J = 1 ^ and J = 2 ^ 1 
transitions are listed for DCO^, while unsplit transition frequencies from J = 2 up to 
J = 14 are reported for all three isotopomers. It is worth mentioning that the quoted 
uncertainties in Table El are la errors derived from the least-squares variance-covariance 
matrix for the fitted spectroscopic parameters. 

As for the accuracy of the measured transition frequencies, the sources of systematic 
error are of two types: those due to the experimental set-up and procedures, and those 
related to the physics of molecules in the experimental environment. 

The sources of the first kind include: accuracy of the frequency standard, harmonic dis- 
tortion of the frequency modulation, integration due to slow lock-in time constant, asym- 
metry in the line shape due to etalon effect in the absorption cell. The well tested per- 
formances of the spectrometer l|Cazzoli fc Dore 1990al |Puzzarim7 Cazzoh fc Dore 2002| ) 
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Table 2. Rotational transition frequencies and spectroscopic constants of DCO 



J'(F') 


J(F) 


observed/MHz 


obs.-calc./kHz 


uncertainties /kHz'' 


1(0) 


0(1) 


72 039.2413 


-0.1 


7.3 


1(2) 


0(1) 


72 039.3028 


-0.3 


6.1 


1(1) 


0(1) 


72 039.3504 


-0.2 


6.7 


2 


1 


144 077.2890*" 


3.9 


5 


3 


2 


216112.5822 


2.1 


5 


4 


3 


288 143.8583 


0.5 


5 


5 


4 


360 169.7783 


-1.2 


5 


6 


5 


432 189.0052 


-1.3 


5 


7 


6 


504 200.1999 


-0.5 


5 


8 


7 


576 202.0239 


1.1 


5 


9 


8 


648193.1357 


0.0 


5 


10 


9 


720172.2024 


1.1 


10 


11 


10 


792 137.8811 


-1.2 


10 




z 1.3 








Constant'* 


this work 


previous" 


DC*^N 


B/MHz 


rotational 


36 019.76765(14) 


36 019.784(2) 




Dj/kHz 


quartic c. d/ 


55.7960(22) 


55.71(12) 


54.3924(15)*^ 


Hj/Rz ; 


jcxtic c. d/ 


0.054(11) 




0.0633(12)*^ 


egg/kH 


z quadrupole 


147.8(35) 




200.6(5)'' 


C//kHz 


spin rotation 


-1.59(78) 




-1.9(1)'' 



" Uncertainties estimated as explained in the text 

^ Not included in the fit, see text 

" (E(obs.-calc.)Vde grees of freedom) ' 

"* Standard errors are reported in parentheses in units of the last quoted digits 

" Bogey et aL llMTI 

* centrifugal distortion 

^ IBriinken et al. 20041 

'' ICazzofi fc Dore 1990bl 

and a judicious choice of the lock-in time constant (less than 1/100 of the rising time 
during the line recording) allowed all these error sources, except the last, to be rendered 
negligible: Figure |21 shows, in fact, a pair of spectra of the J = 5 <— 4 transition of 
DCO^with opposite asymmetry. The shift of the apparent center frequency, however, 
may be accounted for by a full profile analysis with a polynomial function describing 
the baseline and a dispersion term included; this procedure gives an accurate line center 
UDore 2003|l . 

The second kind of error in the line center includes: Doppler shift due to ion drift 
in the discharge, gas flow through the cell due to pumping, pressure shift. The first two 
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Table 3. Rotational transition frequencies and spectroscopic constants of D^'^CO 



J' 


J 


observed/MHz 


obs.-calc./kHz 


uncertainties /kHz'' 


2 


1 


141465.1331'' 


3.4 


5 


3 


2 


212 194.4920 


1.9 


5 


4 


3 


282 920.0055 


0.3 


5 


5 


4 


353 640.3923 


-1.3 


5 


6 


5 


424 354.3718 


-1.8 


5 


7 


6 


495 060.6646 


0.4 


5 


8 


7 


565 757.9861 


2.0 


5 


9 


8 


636 445.0518 


-0.8 


5 


10 


9 


707 120.5885 


-0.7 


15 


11 


10 


777 783.3159 


2.3 


50 




z 1.8 








Constant"* 


this work 


previous^ 


Ql3Ql5j^f 


B/MHz 


rotational 


35 366.70968(21) 


35 366.712(11) 




Dj/kHz 


quartic c. d.^ 


53.4083(40) 


53.20(53) 


52.26535(96) 


Hj/Hz ; 


jextic c. d.^ 


0.048(24) 




0.06074(73) 



" Uncertainties estimated as explained in the text 

^ Not included in the fit, see text 

" (E(obs.-calc.)Vde grees of freedom) ' 

Standard errors are reported in parentheses in units of the last quoted digits 
" Bogey et al. HMH 
^ IBriinken et al. 20041 
^ centrifugal distortion 

causes of line shift should be negligible in the present experimental conditions, and are, 
anyway, largely suppressed by the double-pass arrangement. As for the frequency shift 
by Ar pressure (about 4 mTorr in this experiment), Buffa et al. H1994|l have shown in 
the case of HCO^that it is most significant for the J = 1 ^ transition; an estimate for 
DCO^broadened by Ar (IBiiffa 20n4|l indicates that, among the rotational lines considered 
here, only the J = 2 <— 1 center frequency should be slightly shifted (about 3 kHz); 
therefore it has been excluded from the fit. 

In the present case, there are also the hyperfine contributions of the D and ^^C nuclei 
which may affect the accuracy of the determined transition frequencies. These effects 
are more significant for the lower J transitions, therefore the most affected among those 
included in the fit should be the J = 3 ^ 2. Its hf structure was predicted assuming 
the hf constants from the present work for D and from Schmid-Burgk et al. H2004|l 
for ^^C; then, to simulate the blended experimental profile, modulated Voigt functions 
(see IDore 2003| were summed at each hyperfine frequency with the relative intensity 
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Table 4. Rotational transition frequencies and spectroscopic constants of DC^^O 



J' J 


observed/MHz 


obs.-calc./kHz 


uncertainties /kHz'' 


2 1 


137 653.5239'' 


3.1 


10 


3 2 


206 477.2402 


-0.9 


5 


4 3 


275 297.3150 


1.5 


5 


5 4 


344112.5213 


-1.1 


5 


6 5 


412 921.6530 


0.4 


5 


7 6 


481 723.4893 


0.3 


5 


8 7 


550 516.8156 


-1.7 


10 


9 8 


619 300.4281 


4.4 


35 


10 9 


688 073.0967 


1.9 


45 


a},JkUz 2.5 








Constant'* 


this work 


previous^ 




B/MHz rotational 


34 413.78556(18) 


34 413.798(1) 




Dj/kRz quartic c. d/ 


50.6704(44) 


50.29(4) 




Hj/YLz sextic c. d.* 


0.109(34) 







^ Uncertainties estimated as explained in the text 

^ Not included in the fit, see text 

" (E(obs.-calc.)Vde grees of freedom) ' 

Standard errors are reported in parentheses in units of the last quoted digits 
" Bogey et aL llMTI 
' centrifugal distortion 

as weighting factor; finally, the center frequency of the synthetic profile was measured 
to check for a shift with respect to the assumed unperturbed J — 3 2 transition 
frequency. From this procedure carried out for all three isotopomers, it turned out that 
the apparent line center is shifted at most by 0.3 kHz from the unperturbed center 
frequency: this indicates that the hyperfine structure is irrelevant in this case, and a 
fortiori for the higher J transitions, as seen in Sect. 

4. Conclusions 

The present paper has shown how radioastronomical observations together with labora- 
tory work are needed to determine with high precision the spectroscopic parameters of 
molecular species. In the laboratory, it has been possible to measure the frequencies of 
DCO+, Di^CO^, and DCi^O^ lines up to 792 GHz with a fairly high accuracy. New 
values of the rotational, quartic, and sextic distortion constants have been determined 
for the three isotopomers with such a precision to allow to predict spectra up to 1000 
GHz with an accuracy of at least 1 part in 10^ for DCO^, 5 parts in 10® for D^^CO^, 
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Fig. 3. Spectra of the J = 4 <— 3 transition of DCO+ with opposite hne asymmetry 
due to different standing wave pattern in the ccU: the 1601 points of each spectrum were 
recorded for a total time of 39 s at a rate of 0.8 MHz s^^ with a lock-in time constant 
of 1 ms. Their profile has been fit to a model Galatry profile with the inclusion of a 
dispersion term to recover the line frequency. 



and of 1 part in 10^ for DC^^O^. Previous frequency estimates were much less accurate 
and did not allow to make accurate kinematic studies in dense molecular cloud cores. 

Using the IRAAI 30ni antenna, we detected for the first time the hyperfine structure 
of the DC0^(1 0) line toward the quiescent starless core L1512, in Taurus. This is the 
first time that the hyperfine structure due to a dcuteron has ever been observed; this 
gave us the possibility of obtaining the hyperfine parameters of DCC^. The DCO^(l-O) 
frequency has also been estimated comparing the observed line with a high sensitivity 
spectrum of N2H+(l-0), and using the new value of the N2H+(l-0) line as calculated by 
Dore et al. (2004). This value is however still uncertain by 6 kHz because of astronomical 
errors, and new sub-Doppler laboratory work on N2H+ is planned to improve this preci- 
sion. The estimated DCO'''(1-0) hf frequencies were included with their uncertainties in 
a weighted least-squares fit of the laboratory data to derive accurate spectroscopic con- 
stants. Prom the fitted constants, the hyperfines have been reestimated with an accuracy 
of better than 2.5 parts in 10^. 
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Table 5. Recommended line frequencies'^ (MHz) of DCO+, 0^300+, and 001^0+ 



J'(F') 


J{F) 


DCO+ 


D"CO+ 




1(0) 


0(1) 


72 039.2414(18) 






1(2) 


0(1) 


72 039.3031(9) 






1(1) 


0(1) 


72 039.3506(14) 






2(1) 


1(1) 


144077.2144(18) 






2(1) 


1(2) 


144077.2619(29) 






2(3) 


1(2) 


144077.2804(10) 






2(2) 


1(1) 


144 077.2851(5) 






2(1) 


1(0) 


144 077.3237(15) 






2(2) 


1(2) 


144 077.3326(23) 






3 


2 


216112.5800(6) 


212 194.4901(9) 


206 477.2411(7) 


4 


3 


288 143.8577(6) 


282 920.0052(9) 


275 297.3135(6) 


5 


4 


360169.7795(6) 


353 640.3936(8) 


344112.5224(6) 


6 


5 


432 189.0066(5) 


424 354.3736(7) 


412 921.6526(6) 


7 


6 


504 200.2004(6) 


495 060.6642(8) 


481 723.4890(8) 


8 


7 


576 202.0228(7) 


565 757.9841(9) 


550 516.8173(13) 


9 


8 


648 193.1356(8) 


636 445.0526(11) 


619 300.4237(31) 


10 


9 


720 172.2013(11) 


707 120.5892(28) 


688 073.0947(72) 


11 


10 


792 137.8823(22) 


777783.3135(67) 


756 833.618(14) 


12 


11 


864 088.8419(47) 


848 431.946(13) 


825 580.782(26) 


13 


12 


936 023.7433(88) 


919 065.206(24) 


894 313.375(43) 


14 


13 


1007 941.251(15) 


989 681.816(40) 


963 030.188(68) 


15 


14 


1 079 840.028(24) 


1 060 280.497(62) 


1031730.011(103) 



The la uncertainties are reported in parentheses in units of the last quoted digits 



The hyperfine structure of the DCO~''(1-0) line needs to be taken into account in the 
analysis of astronomical data, to avoid overestimates of the line width by more than a 
factor of 2 and for a correct interpretation of the line profile. 
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